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A METHOD FOR THE CALCULATION OF LARGE NUMBERS OF 
DIPOLE AND QUADRUPOLE TRANSITION PROBABILITIES 

By Lawrence P. Shomo, Goetz K. Oertel,* 
and Charles S. Freert 
Langley Research Center 

SUMMARY 

A computer program is presented which selects allowed transitions and calculates 
dipole and quadrupole transition probabilities for transitions with LS coupling and no 
equivalent electrons, based on an extension of the Coulomb approximation formalism to 
quadrupole and higher multipole transitions. Absorption oscillator strengths or f-values 
calculated by (1) the self-consistent-field method, (2) the scaled Thomas-Fermi method, 
(3) the Coulomb approximation method, (4) the variational method, and (5) the effective 
charge method for singlet and triplet transitions in neutral helium are presented and 
compared. The Coulomb approximation f-values calculated with the present computer 
program are found to be in good agreement with the results obtained by the more sophisti- 
cated methods. 


INTRODUCTION 

The present status of the knowledge of transition probabilities or of absorption 
oscillator strengths (f-values) is poor compared with that of the wavelengths of atomic 
spectral lines. The f-values are known only for a small number of the lines for which 
wavelengths are known, and the accuracy of the available values is usually poor. 

The calculation of f-values requires a knowledge of the wave functions of the two 
states involved in the transition. It is necessary to employ approximations to determine 
the wave functions for nonhydrogenic atoms or ions since Schrodinger's equation can be 
solved analytically only for the simple one-electron atom case. Several techniques that 
provide approximate wave functions have been developed. The most sophisticated tech- 
nique of these is the self -consistent-field method (refs. 1, 2, and 3). Provided that 
exchange and correlation effects are considered, f-values can be obtained which should 
not deviate from the true values by more than 20 percent. Solution by numerical 
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techniques on a computer is almost mandatory because of the complexity and iterative 
nature of the self-consistent -field calculations. The computing time still becomes 
excessive when very complicated systems (many-electron systems) are considered. 

The scaled Thomas-Fermi method (ref. 4) provides a useful alternative. This 
method is not subject to some of the limitations of the Coulomb approximation method 
and is thus more generally applicable. Its accuracy is, of course, limited by the validity 
of the Thomas-Fermi model of the atom's or ion's core (inner electrons plus nucleus). 
The adjustment of a scaling factor, however, to provide energy values consistent with the 
measured values allows for some compensation for the model's inadequacies. Even for 
relatively simple systems, the computation time for the scaled Thomas-Fermi method is 
less than that for the self-consistent-field method by at least one order of magnitude. 

The Coulomb approximation method (ref. 5) is applicable to many lines of interest, 
and for those transitions where its assumptions are met, it gives f-values which are in 
close agreement with the self-consistent-field results. Even for relatively simple sys- 
tems where the time required for self-consistent-field calculations is still relatively 
short, the Coulomb approximation method will result in a reduction of computing time by 
about three orders of magnitude. 

Radial matrix elements for the Coulomb approximation method can be fully charac- 
terized by the effective principal quantum numbers and the individual electron orbital 
angular momentum quantum numbers of the two levels involved in the transition. Thus, 
generalized tables can be assembled in terms of these quantum numbers alone. In this 
sense the Coulomb approximation method has a distinct advantage over the other two 
methods. 

The Coulomb approximation method is attractive because of its simplicity and wide 
range of applicability. Oertel and Shomo (ref. 6) extended the Coulomb approximation 
formalism to the multipole case and presented extensive tables for computing Coulomb 
approximation radial matrix elements for dipole transitions s-p, p-d, d-f, f-g, g-h, and 
h-i, and for quadrupole transitions p-p, d-d, f-f, g-g, h-h, i-i, s-d, p-f, d-g, f-h, g-i, and 
h-j. The extension of the Coulomb approximation method to higher multipole transitions 
was deemed desirable because of the need for knowledge of forbidden electric-multipole 
transition probabilities in astrophysics and because large numbers of matrix elements, 
both dipole and quadrupole, are required for computations of the Stark broadening of 
isolated ion lines (ref. 7). 

The purpose of this paper is to present a FORTRAN IV computer program which 
selects allowed transitions and calculates Coulomb approximation line strengths and tran- 
sition probabilities for dipole and quadrupole transitions with LS coupling and no equiv- 
alent electrons, based upon the formulae presented in reference 6, and to compare the 
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Coulomb approximation f-values with results obtained from more sophisticated calcula- 
tions for neutral helium. A listing of the computer program is presented in appendix A. 
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SYMBOLS 

dipole or quadrupole transition probability, second - 1 
coefficients generated from recurrence relation (eq. (8)) 

Bohr radius, centimeters 

velocity of light in vacuum, centimeters/second 

measured binding energy of orbit with quantum numbers n* and l, 
centimeter - ! 

elementary charge 

absorption oscillator strength, dimensionless 

statistical weight of level a and level a ' , respectively, dimensionless 
Planck's constant/27r, erg-seconds 
radial factor, (-l)^ 5 ” I s 
radial transition integral 

total atomic angular momentum quantum number, dimensionless 
summation indices, dimensionless 

total atomic orbital angular momentum quantum number, dimensionless 

change in total atomic orbital angular momentum quantum number for a 
transition 

total atomic orbital angular momentum quantum number of parent configura- 
tion, dimensionless 
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l,V individual electron orbital angular momentum quantum numbers, 

dimensionless 

Z> larger of the individual electron orbital angular momentum quantum numbers 

of two states involved in transition, dimensionless 

Z< smaller of the individual electron orbital angular momentum quantum numbers 

of two states involved in transition, dimensionless 

n,n' effective (noninteger) principal quantum numbers, dimensionless 

n* integer principal quantum number, dimensionless 

An* absolute magnitude of change of principal quantum number in a transition, 

dimensionless 

R(n,Z),R(n',Z') normalized radial wave functions, atomic units 

^line line factor, dimensionless 

Rmult multi plet factor, dimensionless 

Roo Rydberg constant, centimeter - 1 

r magnitude of position vector, atomic units (aq) 

S total atomic spin quantum number, dimensionless 

S(line) relative line strength, dimensionless 

S(mult) relative multiplet strength, dimensionless 

S(q',Q!') line strength, atomic units ^a^e^ for s = 1, a^e^ for s = 2^ 
s index (1 for dipole, 2 for quadrupole) 

T Q! ,T Q! t term value of level a and level a', respectively, centimeter - l 
Too term value of series limit, centimeter - 1 
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W(abcd;ef) 


Racah coefficient, dimensionless 


Z charge of nucleus minus charge of inner electrons, dimensionless 

Z* effective charge, dimensionless 

a lower level 

a' excited level 

r gamma function 

X transition wavelength, angstroms 

a dipole radial factor 


DEFINITIONS 


The dipole absorption oscillator strength or f-value f(a',a'') 
probability A(a,a!') for a transition between an excited level a' 
can be expressed in terms of the dipole line strength S(a,a') by 


and the transition 
and a lower level a 


and 


1 1 l T o> ~ Ta'l S^a') 

3 gQ, e 2 a2 


(1) 


A (a, a’) 


4 3 

l/e 2 \ c 1 /| T 0! ~ T ot'| ^ S^a*) 

6\Rc/ a 0 I^\ Rco ) e 2 a 2 


( 2 ) 


where aQ is the Bohr radius, T a , T a », g a , and g a » are the term values and statis- 
tical weights of the levels a and at', respectively, Roo is the Rydberg constant, and c, 
e, and h have their usual meanings. 

The line strength can be factored into angular and radial factors, that is, 


S(«,q') = S(q’,«) = S(line)S(mult)cr 2 


( 3 ) 


where the angular part is composed of the relative line strength S(line) and the relative 
multiplet strength S(mult). The dipole radial factor cr2 is related to the radial 


5 



( 4 ) 


transition integral I s (s = 1 for the dipole case) by 



where 


p°° 

I s = \ R(n,Z)R(n’,Z')r s dr (5) 

J 0 

l y is the larger of the individual electron orbital angular momentum quantum numbers l 
and V, and R(n,Z) and R (n',r) are the normalized radial wave functions of the active 
electron in atomic units. Relations for the angular factors of the dipole line strength are 
given in appendix B. 

Analogous expressions for the transition probability and line strength for quadrupole 
transitions are given in appendix C, and they are discussed in detail in reference 6. 

COULOMB APPROXIMATION 


Bates and Damgaard (ref. 5) assert that for most transitions the potential approxi- 
mates closely its asymptotic (r — °°) Coulomb form in the region of dominant contribu- 
tion to the transition integral I s . This is due to the emphasis placed on long radial dis- 
tances by the factor r s in Is, an emphasis which becomes greater with an increase in 
s, that is, higher multipole transitions. At sufficiently large electron-nucleus separa- 
tions, the active electron is outside the core of inner electrons and is acted upon by a 
potential which is very nearly Coulombic. The Coulomb approximation wave functions are 
obtained by solving the standard central-field radial Schrodinger equation using the 
Coulomb potential, subject to the condition that the solutions approach zero at large r. 
The additional constraint that the energy parameter in the Schrodinger equation be identi- 
fied with the experimentally determined binding energy of the level is imposed. Since the 
energy parameter defined in this manner is not, in general, an eigenvalue of the 
Schrodinger equation, the solutions diverge at the origin. 

In terms of the Coulomb approximation wave functions R, the radial transition 
integral I s (s = 1 for dipole, s = 2 for quadrupole, etc.) for transitions between levels 
(n,Z) and (n’,Z') is obtained from 


I s = C R(n,Z)R(n',r)r s dr 
J 0 

nn+n'rrn+n'+ir V C 

= ^ 5 r(n + i + l)r(n - I)r(n' + 1 ' + I)r(n' - i’) > > a k a k , \ 

n n+ V) n ' +lL J £)kM) L JC 


e -rZ(n+n t /nn f ) r n+n , +s-k' 


k 'drj 


( 6 ) 
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where r denotes the gamma function. The principal quantum numbers n and n' are 
noninteger (except for hydrogenic configurations) "effective" principal quantum numbers 
and may be obtained from 


n = z (? 


109737\ 1//2 

T) 


(7) 


in terms of the series limit Too (in units cm-1), the measured term value T, and Z 
the charge of the nucleus minus that of the inner (core) electrons (i.e., Z = 1 for neu- 
tral atoms). The normalization factor proposed by Hartree (ref. 8) for hydrogenic 
radial wave functions was used although it is not accurate for noninteger principal quan- 
tum numbers. The error thus introduced, however, was shown to be quite small for most 
cases of interest (ref. 6). The coefficients a^ were generated from the recurrence 
relation 


na, 


a k = 


k- 

2kZ 


IE 


l(l + 1) - (n - k)(n - k + 1) 


■) 


( 8 ) 


with a.Q = 1. The coefficients a^/ are given by equation (8) with n and l replaced 
by n' and V , respectively, and with apt = 1. The integral over the exponential func- 
tion in equation (6) can be evaluated as 




nn' 


-,n+n’+s+l-k-k' 


Z(n + n') 


r(n + n' + s + 1 - k - k') 


(9) 


In order to avoid the divergence of the wave functions at the origin, the summation is cut 
off by the condition k + k’ <n + n' + s - 2, which amounts to neglecting powers of r less 
than 2 in equation (6). Thus, the basic assumptions of the Coulomb approximation of Bates 
and Damgaard are violated when the contribution to the transition integral is large in the 
core and small at long distances (i.e., when the active electron is in an inner orbital with 
very small amplitude outside the core, or due to cancellation outside the core region). 


COMPUTER PROGRAM 


A FORTRAN IV computer program was written to calculate dipole f-values, and both 
dipole and quadrupole wavelengths, transition probabilities, and line strengths for all of 
the allowed (dipole and quadrupole) transitions for a given element and stage of ionization. 
The program is restricted to levels with LS coupling and no equivalent electrons. It cal- 
culates the radial matrix elements by the Coulomb approximation formalism developed in 
reference 6 and described in the present paper. The program is based upon the equations 
presented in this paper and is listed in appendix A. 


7 



The transition parameters can be calculated for all the allowed transitions in a 
given element and stage of ionization, or the transition parameters can be calculated only 
for certain allowed transitions by restricting the input data to only those levels involved 
in the desired transitions. The program operates as follows. The program selects one 
of the input levels as the initial level and then scans all other levels and applies dipole 
and quadrupole selection rules to determine the allowed final levels. As each final level 
is determined, the required transition parameters are calculated. Then, another level is 
selected as the initial level and the process is repeated. This procedure continues until 
all the input levels have been chosen as the initial level. 

The Coulomb approximation method described in this paper and in reference 6 is 

utilized in subroutine BATES. The output of this subroutine (RSQUAR) is the value of the 

o 

radial transition integral squared, that is, I s . Racah coefficients (ref. 9) are calculated 
in subroutine RACAH according to the relation 


W(abcd;ef) = (-l) a+b+c+d A(abe)A(bdf)A(afc)A(edc) (-1)% + l)l[(/3 - a - b - e)! 

0 

x (0 - b - d - f)»(j8 - a - f - c):(/3 - e - d - c):(a + b + c + d - j8)! 
x(a + d + e+ f- /3).'(b + c + e + f - /3)fj 


where 


A(ijk) = 


(i + j - k):(j + k - i):(k + i - j) ! 
(i + j + k + 1): 


1/2 


( 10 ) 


(H) 


is calculated in subroutine DELTA and where the summation in equation (10) is over those 
integer values of /3 for which the factorials are meaningful. Extensive tables of Racah 
coefficients are available (ref. 10). 

The computer program offers several options controlled by a set of option variables 
as follows. When OPQUA = 1, only dipole transitions will be calculated; but when 
OPQUA = 2, both dipole and quadrupole transitions will be calculated. When OPFVA = 1, 
dipole f-values and dipole transition probabilities will be calculated; when OPFVA = 2, no 
f-values or transition probabilities will be calculated; and when OPFVA = 3, dipole 
f-values and dipole plus quadrupole transition probabilities will be calculated. The line 
strength, wavelength, and the angular frequency are calculated regardless of the value of 
OPFVA. The option controlled by NPROP was introduced to provide an easy means of 
checking the program results by hand calculation. This option provides the entry values 
for the tables of Coulomb approximation radial matrix elements that were presented in 
reference 6. The correct value of the radial matrix element for use in the hand calcula- 
tion can thus be obtained from these tables. When NPROP = 1, these table-entry values 
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will be printed out for each transition; but when NPROP = 2, they will not. The value of 
the option variable NLITO is the upper limit of An* (the absolute magnitude of the dif- 
ference between the integer principal quantum numbers for the two levels involved in the 
transition). The program will not consider any transitions with An* greater than 
NLITO. An upper limit is set on NLITO by statements in the program which require 
NLITO = 9.0. The first input card contains the values of these four option variables 
OPQUA, OPFVA, NPROP, NLITO as required by FORMAT statement 140. 

The second and third input cards contain information defining the atomic system 
under consideration. The second card contains the following input variables: 


ELEMNO atomic number of element 


Z core charge of atomic system 

NOLEV number of input terms (NOLEV is independent of J-values) 

Any number =90 of no-equivalent electron and LS coupling terms can be 
used. 

NOSER number of series in atomic system (the limit NOSER < 10 is contained 

in the program) 


The second input card contains the values of the variables ELEMNO, Z, NOLEV, and 
NOSER as required by FORMAT statement 3. 

There will be NOSER input cards like the third card, that is, one card is required 
for each series. Each card will contain the following input variables: 

CAPSS(I) total atomic spin quantum number of Ith series designation 

CAPLS(I) total atomic orbital angular momentum quantum number of Ith series 

designation 

SLIM(I) series limit (in units cm-1) of Ith series 

Thus, the next NOSER input cards (starting with the third input card) each contain the val- 
ues of the variables CAPSS(I), CAPLS(I), and SLIM(I) for one of the NOSER series as 
required by FORMAT statement 7. 

The next series of input cards define the terms of the atomic system. A sequence 
of two input cards is required for each term. There are NOLEV of these two-card 
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sequences for the atomic system. The first card of the sequence contains the following 
input variables : 


NLITT(I) 

LLITT(I) 

CAPS(I) 

CAPL(I) 

PARITY (I) 

CAPSP(I) 

CAPLP(I) 


principal quantum number of Ith term 

individual electron orbital angular momentum quantum number of Ith term 

total atomic spin quantum number of Ith term designation 

total atomic orbital angular momentum quantum number of Ith term 
designation 

parity of Ith term (PARITY(I) = 0,2 means odd or even parity, 
respectively) 

total atomic spin quantum number of Ith term's parent configuration 

total atomic orbital angular momentum quantum number of Ith term's 
parent configuration 


Thus, the first card of each two-card sequence defining a term contains the values of the 
variables NLITT(I), LLITT(I), CAPS(I), CAPL(I), PARITY(I), CAPSP(I), and CAPLP(I) as 
required by FORMAT statement 6 . 

There are NOJ J-values (total atomic angular momentum quantum numbers) cor- 
responding to a given term, where NOJ is an integer =1 calculated in the program. The 
second card in the two-card sequence contains the energy level values of the term cor- 
responding to each of the term's J-values (i.e., TJ(N,I), N = l,NOJ (of the Ith term)). 
The term values must be listed on the card in order of increasing J-value. TJ(N,I) is the 
energy level (in units cm~l) of the Ith term and Nth J-value. Thus, the second card of 
each two-card sequence defining a term contains the values of the variables TJ(1,I), 
TJ(2,I),. . ., TJ(NOJ,I) as required by FORMAT statement 4. 


In order to clarify the foregoing instructions, a sample input deck is listed at the 
end of the program in appendix A. This sample input deck contains six levels of neutral 
helium, namely, 2 s( 2 S) 3 Sj, 2p(2s) 3 p P j 2 ,3p( 2 S) 3 Pq j 2 , Sd^S^Dj 2 3 , 4p ^ 2s ^ 3p 0,l,2’ 
4d( 2 S) 3 Dj 2 3 (where the superscript 0 indicates odd parity). The term values used are 

those of C. E. Moore (ref. 11). Comment cards have been inserted in the sample input 
deck to explain each input data card. 
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A sample computer program printout follows the sample input deck listing in 
appendix A. This printout is a partial printout for the sample input deck (partial in that 
only the printout is given for the case in which the 2 s(2S)3si level is selected as the ini- 
tial level by the program). The designation INITIAL LEVEL used in the program print- 
out refers to the lower state involved in the transition. The printout for a given initial 
level consists of the initial level identification with appropriate parameters, followed by 
a list of the allowed final levels. This list of final levels consists of (1) the final level 
identification and (2) the calculated transition parameters for that line (provided it is an 
emission line, i.e., X > 0). Further, the final levels in the list are grouped as quadru- 
pole transition levels (with subgroups |AL| = 0, 1, 2 ) or dipole transition levels. 

The output variables are presented with the printout variable name listed first and 
the program variable name (used internally in the program) listed next to it in brackets. 
The definition of the variable is given except when the variable has already been defined 
as an input variable. 

A [aAAPn) dipole or quadrupole transition probability 

ABS. OSCILLATOR STRENGTH (fAAPN] dipole absorption oscillator strength or f-value 
CAPITAL L [CAPL(lj] 

CAPITAL S (CAPS (l[] 

DIFFERENCE [DIFF] - n t 

ELEMENT NUMBER (ELEMNO] 

J [J(N,I]1 Nth J-value of the Ith term 

LAMDA [xlAM] dipole or quadrupole transition wavelength in angstroms 
LITTLE L (LLITT(I)] 

LITTLE N [NLITT(I)] 

N STAR OF LARGER LITTLE L VALUE (STARN(LGr[] n t> (where n is effective 

principal quantum number) 

% 


N STAR OF SMALLER LITTLE L VALUE (STARN(LSM)] 
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PARENT TERM CAPITAL L [CAPLP(I)] 

PARENT TERM CAPITAL S [CAPSP(I)] 

PARITY [PARITY (I)] 

RAA |RAA(I)] (r^ for a given level 

S [Is] - quadrupole 

a o e 

S [SAAPN] dipole 

a 0 e 

WAADP (WAADP) quadrupole transition angular frequency (w = 2irt where f is 

frequency) 

WAAP jWAAP] dipole transition angular frequency (w = 27 rf where f is frequency) 

z[z] 


This computer program was run on a Control Data 6600 computer system in which 
14-place precision is standard. It is recommended that more than 8-place precision be 
employed if possible, since 8 -place precision will give accurate Coulomb approximation 
results only for transitions with n + n' less than approximately 10 (ref. 6). 

Numerous comment cards have been incorporated into the listing of the computer 
program given in appendix A to assist the reader in understanding the functioning of the 
program. 

The sample input deck in appendix A can be used as a test deck by comparing the 
printout generated for the initial level 2 s( 2 S)^Si with the sample printout also given in 
appendix A. 


RESULTS AND DISCUSSION 

The results of calculating f-values for several transitions in neutral helium by 
several different methods are presented in table I. The f-values are averaged over the 
J-values of the initial levels by using the appropriate statistical weights and summed 
over the J-values of the final levels. 
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TABLE I.- ABSORPTION OSCILLATOR STRENGTHS FOR NEUTRAL HELIUM 


Transition 

Wiese, Smith, 
Glennon 
(ref. 13) 

Coulomb 

approximation 

(a) 

Self-consistent-field 

method 

Scaled Thomas 
method 

-Fermi 

Variational 
calculation of 
Goldberg 
(ref. 17) 

Variational 
calculation of 
Hylleraas 
(ref. 18) 

Sum rule 
modified 
calculation 
(ref. 19) 

Coulomb 

approximation 

(b) 

Effective 

charge 

method 

l l S - 2lp 

0.2762 

SP 


'0.2531 1 


0.2719 





0.3555 




3*P 

.0734 

SP 


.06825 


.07203 





.0722 




4 lp 

.0302 

W 


.02795 

> 

.02738 





.0282 




5*P 

.0153 

LS 


.01414 











6*P 

.00848 

LS 


(^.008134, 


.00793 





.0082 




2 J S - 2*P 

.3764 

SP 

.3723 

.3578 


0.342 ' 


0.389 

.392 

0,370 

0.3719 

0.3861 

3*P 

.1514 

SP 

.1492 

.1646 



.242 


.157 

.150 

.156 

.1431 

.2948 

4 lp 

.0507 

W 

.04865 

.0508 



.0684 


.057 

.062 

.0503 

.04970 

.08728 

5*P 

.0221 

DK 

.02218 




.0296 


.0252 

.025 

.0221 

.01936 

.03818 

6*P 

.0128 

T 

.01209 

.01257 


. 

.0157 


.0136 

.012 

.0124 



2*P - 3 X S 

.0480 

W 

.04806 

.04640 


.0278 

1 



.04797 


.04859 

.02304 

4*s 

.00824 

T 

.008579 

.00824 


.00549 



.00686 


.008308 

.004150 

5 IS 

.00308 

CA 

.003213 



.00211 





.003041 

.001519 

6 IS 

.00153 

T 

.001593 

.00151 


,.00105_ 





.001491 

.0007502 

2 l P - 3*D 

.711 

W 

.7104 

.7253 


'.575 


.755 

.7357 


.7006 

.6947 

4*D 

.122 

T 

.1207 

.1205 


.114 


.118 

.132 


.1168 

.1199 

5*D 

.0436 

T 

.04346 

.04308 


.0436 


.0416 

.04313 


.04106 

.04355 

6*D 

.0213 

CA 

.02105 



,.0218. 


.0199 



.02054 

.02115 

2 3 S - 2 3 P 

.5391 

SP 

.5457 

.6002 


.538 


.542 

.559 

.535 

.5319 

.6686 

3 3 P 

.06446 

SP 

.05982 

.05705 


.0638 


.0826 

.052 

.0768 

.06478 

.1876 

4 3 P 

.0231 

W 

.02370 

.02083 


.0256 


.0270 

.031 

.0232 

.02581 

,06502 

5 3 P 

.0114 

DK 

.01143 



.0124 


.0123 

.013 

.0114 

.01022 

.02999 

6 3 P 

.00608 

DK 

.006372 



.00694 

.00665 

.007 

.0061 



2 3 P - 3 3 S 

.0693 

W 

.06813 

.07639 


.0692 



.07662 


.07204 

.02700 

4 3 S 

.0118 

T 

.01014 

.01159 


.0105 



.00520 


.008939 

.004288 

5 3 S 

.00365 

CA 

.003581 



.00376 




.003159 

.001488 

6 3 S 

.00176 

CA 

.001720 



.00182 




.001527 

.0007089 

2 3 P - 3 3 D 

.609 

W 

.6168 

.6234 


.611 


.553 

.6530 


.6089 

.6969 

4 3 D 

.125 

T 

.1244 

.1229 


.123 


.129 

.0792 


.1192 

.1325 

5 3 D 

.0474 

T 

.04759 

.04668 


.0472 


.0512 

.04951 


.0455 

.04958 

6 3 D 

.0215 

CA 

.02374 



.0236 


.0260 




.02446 


Calculated by computer program in appendix A. 

Calculated by linear interpolation in original tables of Bates and Damgaard (ref. 5). 


Notation: 

SP Results of variational calculations from reference 14. 

W Results of variational calculations of A. W. Weiss presented in reference 13. 
LS Results of Low and Stewart quoted in reference 15. 

CA Coulomb approximation results presented in reference 13. 

DK Results of reference 19 modified by Wiese, Smith, and Glennon (ref. 13). 

T Results of reference 16 modified by Wiese, Smith, and Glennon (ref. 13). 

£ J Stewart-Rotenberg scaling factor is less than zero (see appendix D). 

£ y Coefficient of fractional parentage used. 
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Perhaps the two most extensive sources of information on experimental and cal- 
culated f-values available to date are two National Bureau of Standards publications: 
namely, an extensive bibliography on atomic transition probabilities by Glennon and Wiese 
(ref. 12) and tables of atomic transition probabilities for hydrogen through neon selected 
after extensive evaluation of all available material by Wiese, Smith, and Glennon (ref. 13). 
The f-values chosen in reference 13 as the best values available are presented in the sec- 
ond column of table I. The results in this column include (1) those of A. W. Weiss calcu- 
lated by a variational method in a manner similar to that used by Schiff and Pekeris 
(ref. 14), (2) Coulomb approximation results calculated by Wiese, Smith, and Glennon with 
a computer program made available to them by H. R. Griem, and (3) those of Low and 
Stewart (unpublished but quoted by Dalgarno and Stewart in ref. 15). Coulomb approxi- 
mation results calculated by the computer program listed in appendix A, which uses the 
Coulomb approximation relations given in the present paper, are presented in the third 
column of table I. The self-consistent-field results including exchange effects (Hartree- 
Fock method) are those of Trefftz and associates (ref. 16). In addition to exchange 
effects, they attempt to take into account certain correlation effects, noncentral nature of 
the real system, and configuration mixing through the introduction of parameters of inter- 
electron distance. The scaled Thomas-Fermi results were made available by John 
Cooper of the Joint Institute for Laboratory Astrophysics, Boulder, Colorado. These 
results were calculated by the scaled Thomas-Fermi method of Stewart and Rotenberg 
(ref. 4). (See appendix D for discussion.) Also presented are the f-values obtained by 
the variational calculations of Goldberg (ref. 17) and Hylleraas (ref. 18) and f-values 
modified by Dalgarno and Kingston (ref. 19) such that they satisfy four different sum 
rules. Coulomb approximation results calculated by linear interpolation in the original 
Coulomb approximation tables of Bates and Damgaard (ref. 5) are listed in the ninth col- 
umn of table I. Finally, results obtained by an effective charge method are given. The 
effective charge method utilizes eigenfunctions of the Schrodinger equation for hydrogen 
with an effective charge Z*, where Z* is chosen such that the hydrogenic energy eigen- 
value equals the measured energy value E n * ^ of the system, that is, 


Z 


*2 


Roc 


(12) 


where n* is the integer principal quantum number. 

A problem encountered in attempting to compare f-values calculated by several 
methods is the lack of knowledge of the true f-value and the difficulty in estimating the 
error introduced by the different approximations. In making comparisons among the 
present tabulated values, the values presented in reference 13 may be used as a general 
standard, although in certain cases, which are discussed, some of the present values may 
be more accurate. 
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For singlet transitions, the scaled Thomas-Fermi method is not rigorously appli- 
cable. The comparison of the new Coulomb approximation results (third column of 
table I) with the values from reference 13 is quite favorable. The good agreement of the 
Coulomb approximation f-values for triplet transitions with the values from reference 13 
is immediately apparent. For triplet states, the scaled Thomas-Fermi method is appli- 
cable, and when significant differences exist between the scaled Thomas-Fermi results 
and the Coulomb approximation results, one can expect that these differences are due to 
significant contributions to the transition integral by the core region of the wave functions 
involved. As the levels involved in the transition become more highly excited levels, one 
expects the difference between the scaled Thomas-Fermi and Coulomb approximation 
results to decrease due to the decreasing contribution of the core regions for such wave 
functions. The results for the 2^p - n^D series in table I clearly show this trend with 
the 2^P - 6^D transition having the smallest difference for the series (about 0.6 percent 
compared with a 1.1 percent maximum for the series). However, this trend is not evident 
for the 2^S - n^P and 2^P - n^S series where the orbits involved are generally closer to 
the core. The average difference for the 2^S - n^P series (about 5.9 percent) is larger 
than that of the 2^p - n^S series (about 3.9 percent) as was to be expected due to the 
presence of the 2^S level with large core penetration in each transition integral of the 
former series. For triplet transitions with a significant difference between the results 
of the two methods, the scaled Thomas-Fermi results can be expected to be the more 
accurate. The fact that some of these Coulomb approximation results are in better 
agreement with the values from reference 13 may not be significant. 

The effective charge results are generally poor. This failure of the effective 
charge method can be attributed to the fact that the method fits a hydrogenic wave function 
for all values of r, including the generally insignificant core region. Although the Bates- 
Damgaard Coulomb approximation is poor in the core region, it is certainly best for 
larger values of r. 

Reference 13 gives Coulomb approximation results as the best available values for 
several of the transitions listed in table I. For these cases, the present Coulomb approx- 
imation results (presented in the third column of table I) should be more accurate since 
the values from reference 13 were obtained by interpolation in the tables of Bates and 
Damgaard. Some feeling for the size of the errors introduced by such interpolation can 
be gained by comparing the present Coulomb approximation results in the third column of 
table I with the results in the ninth column of the table which were calculated by a com- 
puter program that interpolated linearly in the Bates -Damgaard tables. The results 
often differ in the second significant figure. 

It was to be expected that the Coulomb approximation results for the series of tran- 
sitions involving the ground state of neutral helium would, in general, be least accurate 


15 



due to the large contributions to the transition integral from the core region of the (two 
equivalent electron) ground state. The values from reference 13 for this series were all 
calculated from variational wave functions (one partial exception to this statement is 
the results of Low and Stewart (l^S - 5lp, l^S - 6^P; see ref. 15) who used modified 
hydrogenic wave functions for the upper levels and a variationally determined wave func- 
tion for the ground state). The Coulomb approximation results for the 1*S - n*P series 
were obtained by multiplying the computer program results by the proper coefficients of 
fractional parentage, since the computer program does not handle equivalent electron 
transitions in its present form. The differences between the Coulomb approximation 
results and the values from reference 13 are indeed generally larger for this series of 
transitions than for any of the other series presented. The results for the l^S - 2^P 
transition differ by about 8.7 percent, whereas those for the higher transitions in the 
series differ by not more than about 7.9 percent (the percent differences quoted are rela- 
tive to average value of the two results). The self-consistent-field results are in better 
agreement with the values from reference 13 for the two lowest members of the series. 
For the other members of the series, the Coulomb approximation results are in better 
agreement than are the more sophisticated self-consistent-field results. 

CONCLUDING REMARKS 

The agreement of the Coulomb approximation f-values (absorption oscillator 
strengths) with the presently accepted f-values for neutral helium (i.e., the values from 
U.S. Dep. Com. NSRDS-NBS 4) is indeed satisfactory except for those transitions which 
involve the ground state. For most situations, where inaccuracies of about 10 percent in 
f-values can be tolerated, the Coulomb approximation f-values for the l^S - n*P series 
are also adequate. While it should be understood that the Coulomb approximation may 
not work as well for other atoms, its nearly uniform success for neutral helium is 
impressive. The present results suggest that in situations where systematic calcula- 
tions of large numbers of f-values or related quantities are required, the Coulomb 
approximation method should be used whenever applicable and be checked against the 
scaled Thomas-Fermi method when significant core contributions to the transition inte- 
gral are expected. This procedure should reduce computing times considerably. The 
only alternative method for systematic calculations presently available is the self- 
consistent-field method which usually entails long computing times and which probably 
does not result in significant improvements over the scaled Thomas-Fermi method when 
the latter method is applicable (i.e., for those transitions which result in positive values 
for the Stewart and Rotenberg scaling factor). The variational method does not provide a 
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practical alternative for systematic calculations, since it involves time-consuming trial 
and error calculations and becomes exceedingly complex for excited states. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 31, 1970. 
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APPENDIX A 


FORTRAN IV PROGRAM FOR CALCULATING DIPOLE AND 
QUADRUPOLE TRANSITION PROBABILITIES 


A listing of the computer program described in the main text is given in this 
appendix. A sample input case and the resulting output listing are included. 


Computer Program Listing 

PROGRAM TRAPRO ( INPUT,0UTPUT,TAPE5= INPUT, TAPE7=0UTPUTJ 
INTEGER ELEMNO»OPFVA , OP QUA 
REAL NLITT»LLITT,J,NOJS» NLITO 

DIMENSION NLITT (90) , LLI TT( 90 1 ,CAPS(90) ,CAPL(90) ,CAPSS( 101 »SLIM( 10) 
1 ,TLIM<90) ,N0JS<901, J (10,90) ,T J ! 10 ,90 1 , CAPLP ! 90 1 , A VET J < 90 1 , S TARN < 90 
2), PARITY! 90 1 ,CAPSP(90J , CAPL S ( 10 1 , R AA ( 90 I 
DATA RYDBGD, C»HBAR, ALPH A, AOBOHR , P I /l. 097373 1E+05 , 2.997925E+ 10, 1 . 05 
1 450 E-27.7.29720E— 03, 5.29167E-09,3. 1415927/ 


FORMAT STATEMENTS 


3 FORMAT ( 15, F5. 1,2151 

4 F0RMAT(6F12.4) 

6 FORMAT (7F5.11 

7 FORMAT ( 2F5. 1 , FI 5.4 1 

101 FORMAT (1X33HDELT A L EQUALS 0 QUADRUPOLE TERMS//! 

102 FORMAT! 1X9HLITTLE N=,F5.0,11H, LITTLE L=,F5.0,12H, CAPITAL S=,F5.1 
1 , 12H , CAPITAL L=,F5.1,9H, P ARI TY= , F3 . 0 1 

103 FORMAT ( 1X24H IN I TI AL J VALUE J=,F5.1,10H, RAA= , F10 .3// 1 

104 FORMAT 11X72H FINAL J VALUES S WAADP 

1 LAMDA 1 

105 FORMAT ( 4XF5 . 1 , 12XE 1 2.4 , 5XE 12. 4, 5XF16. 1 1 

106 FORMAT (1X33HDELT A L EQUALS 1 QUADRUPOLE TERMS//! 

107 FORMAT ! 1 X8H J VALUE=, F4. 1,5X26H ABSORPTION LINE ****♦***#*//) 

108 FORMAT ( 1X33HDELT A L EQUALS 2 QUADRUPOLE TERMS//! 

109 FORMAT C 1X13HINITIAL LEVEL1 

110 FORMAT (1X11HFINAL LEVEL! 

111 FORMAT ( 1 X29HP ARENT TERM CAPITAL S=,F5.1,12H, CAPITAL L=,F5 

1.11 

112 FORMAT 1 1 X12HD I POLE TERMS// i 

113 FORMAT ( 1X32HN STAR OF LARGER LITTLE L VALUE=,F8.31 

114 FORMAT ( 1X33HN STAR OF SMALLER LITTLE L VALUE=,F8.31 

115 FORMAT ! 1 X60H*****B0 TH STATES HAVE BEEN ASSIGNED THE SAME TERM VALU 

IE*****) 

116 FORMAT ! 1X72H FINAL J VALUES S WAAP 

1 LAMDA 1 

118 FORMAT (1H1////1X15HELEMENT NUMBER, 1 5, 9H, WITH Z=,F4.01 

119 F0RMAT14XF5.1, 12XE12 .4//! 
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125 FORMAT < I XI 1HDI FF ER ENCE=, F8 .3 I 

126 FORMAT ( 1X3HA= ,E12.4//) 

140 F0RMAT(3I10,F10.1) 

511 FORMAT ( 13H NOSER GT 10) 

717 FORMAT (26H ABS. OSCILLATOR STRENGTH=E 12 . 4, 3H A=E12.4//I 

833 FORMAT ( //24H ERROR, WRONG TERM LIMIT) 

834 FORMAT ( / / ///3HEL=I3»4HLEV=I3»3H N=F5. 1 , 3HSL=F5. 1 , 3H S=F 5. 1 , 3HBL=F5 
1.1, 3H J=F5. 1 ) 

835 FORMAT!//// /4H T J=F 1 2. 2 , 10X6H TLIM=F12.2) 

916 FORMAT (Z///5H*****// /// ) 


DATA READ IN AND ORGANIZED 


READ! 5, 140) OPQUA,OPFVA,NPROP,NLITO 
READ! 5 ,3 ) ELEMNO , Z, NOLEV , NOSER 
100 DO 67 1=1,90 
TLIMII )=0. 

IF! I.LT.ll) SLIM! I) = 0. 

67 CONTINUE 
I *0 

IF! NOSER. LE. 10) GO TO 510 

PRINT 511 

STOP 

510 CONTINUE 
1 = 1+1 

READ (5,7) !CAPSS(IltCAPLS( I)»SLIM(I),I=1* NOSER) 

DO 10 1=1, NOLEV 

RE AD (5,6) NLITT! I),LLITT(I) ,C APS! I ) , CAPL! I ) , PARITY (I ) ,CAPSP( I) ,CAP 
1LP! I) 

DO 9 M=l, NOSER 

IF! ABS(CAPSP( I )-CAPSS(M ) ) .LT..5.AND. ABS ( C APL P ( I ) -CAPLS ( M) ) .LT. .5) 

1 TL I M ( I )=SLIM( M) 

9 CONTINUE 

THE NUMBER OF J VALUES OF A TERM CALCULATED 

DIFF = ABS(CAPL(I)-CAPS( I )) 

SUM=CAPL ( I)+CAPS( t) 

NOJS! I )=SUM-DI FF+1. 

NOJ=NO JS ( I l+.OOl 

J ARRAY MUST BE IN ASCENDING ORDER 

READ! 5,4) ( T J ( N , I ) , N = 1 , NOJ ) 

DO 801 N= 1 » NO J 
XN=N 

801 J(N, I)=DIFF+XN-. 99999999 

LOOPS 871,872 GENERATE MISSING TJ VALUES 

SAV ET =0. 

DO 871 NK=1» NOJ 

IF! TJ(NK, D.LT.TLIM! I)) GO TO 871 
PRINT 833 

PRINT 834, NKt I, NLITT! I ) *LL ITT! I ) , CAPSS(I) , CAPLS! I) ,J(NK,I) 

PRINT 835, TJ(NK,I) ,TLIM( I) 
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871 CONTINUE 

DO 872 NK=1 , NO J 

IF(TJ(NK»I1.LT.1.)TJ(NK»I 1=SAVET 

872 CONTINUE 
10 CONTINUE 

C 

C 

C INITIAL LEVEL SELECTED AND ITS PARAMETERS CALCULATED 

C 

C 

DO 12 1=1, NOLEV 

X=0.0 

Y=0 .0 

NOJ=NOJS( I) 

C 

C CALCULATION OF AVERAGE TERM VALUE 

C (THE AVERAGE TERM VALUE IS THE AVERAGE OF THE ENERGY LEVEL VALUES, 

C EACH WEIGHTED BY THE APPROPRIATE STASTICAL WEIGHT, COMPOSING THE 

C TERM. THIS AVERAGE TERM VALUE IS USED TO CALCULATE THE EFFECTIVE 

C PRINCIPAL QUANTUM NUMBER FOR THE TERM. ) 

C 

DO 600 N= 1 » NO J 
X=X+TJ(N»I)*(2.*J(N, II +1.1 
Y=Y+2.*J(N,I 1+1. 

600 CONTINUE 

AVETJ ( I 1 = X/Y 
C 

C EFFECTIVE PRINCIPAL QUANTUM NUMBER CALCULATED FOR EACH TERM 

C ( SEE EQUATION (7) I 

C 

STARN( I l=SQRT< Z**2*RYDBGD/( TLIM( I) -AVETJ ( I 1 ) ) 

CALL BATES(STARN( I) , LLITTII I ,STARN ( 1 1 , LLI XT ( 1 1 ,2. , Z, RAAI I 
RAA(I) - SORT ( RAAI I 
12 CONTINUE 
C 

C INITIAL LEVEL SELECTED 

C 

DO 915 1=1, NOLEV 
NOJ=NOJS ( I ) 

DO 15 N*1 » NOJ 

IF(NLITT( H.GT.9.51 GO TO 15 
C 

C INITIAL LEVEL PRINTED OUT 

C 

PRINT 118, ELEMNO, Z 
PRINT 109 

PRINT 102, NLITT( II , LL I TT ( I ) ,CAPS ( II , CAPL( I), PARITY (I) 

PRINT ill, CAPSP ( I) , CAPLP( I I 
PRINT 103, J ( N, I ) *R A A( I I 
C 
C 

C DELTA CAPITAL L = 0 QUADRUPOLE TERMS 

C 

C 

I F( OPQUA.EQ.il GO TO 760 
PRINT 101 
F=2 • 
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QUADRUPOLE SELECTION RULES APPLIED IN NEXT 8 STATEMENTS 

IF(ABS(PARITY(I)— PARITY (MN)).GT«. 5 )G0 TO 17 
IF(ABS(CAPS(I) -CAPS ( MN) ) .GT • . 5) GO TO 17 
IF( ABSCCAPH I l-CAPL ( MN) ).GT..5) GO TO 17 
I F < ABS(LLITT(I)+LLITT(MN)).LT..5) GO TO 17 
IF(ABS(CAPSP( II —CAP SP(MN)).GT..5) GO TO 17 
IF* ABS(CAPLP(I)-CAPLP(MN)).GT..5) GO TO 17 
IF( ABS(LLITT(I)-LLITT(MN) ). GT.2.5) GOTO 17 

I F( 0.5.LT.ABSCLLITT( I)- LL ITT( MN ) ) . AND. ABS ( LL I TT(I) -LLI TT ( MN )) . LT . 1 
1.5) GO TO 17 

1F( ABS(NLITTd)— NLITT(MNJ).GT.NLITO) GO TO 17 
FINAL OUADRUPOLE LEVEL IDENTIFICATION PRINTED OUT 
PRINT 110 

PRINT 102 » NLITTIMNI ,LL ITT ( MN ) , CAP S < MN ) ,CAPL(MNI ,PARITY(MN) 

PRINT 111, CAPSP ( MN) ,CAPLP( MN) 

DISTINCTION MADE BETWEEN LEVELS WITH SMALLER AND LARGER LL I TT 
754 B IGSML = LL ITT ( I ) 

IF(LLITT< I)-LLITT(MN).LT.O, ) BI GSML=LLI TT ( MN) 

I F( ABS ( BI GSML-LL ITT ( I)) .LT..5) LGR= I 

IF* ABS(BIGSML-LLITT ( MN) ) .LT..5) LGR* MN 
IF* ABS(BIGSML-LLITT( I) ) .LT..5) LSM=MN 
I F( ABS(BIGSML-LLITT ( MN) ) .LT..5) LSM=I 
DIFF=STARN(LSM)-STARN(LGR) 

IF(NPROP.NE.l) GO TO 905 

COULOMB APPROXIMATION QUADRUPOLE RADIAL MATRIX ELEMENT TABLE 
ENTRY POINTS PRINTED OUT 


PRINT 113, STARN(LGR) 

PRINT 114, STARNCLSM) 

PRINT 125, DIFF 

SQUARE OF QUADRUPOLE RADIAL TRANSITION INTEGRAL CALAULATED 

905 CALL BATEStSTARN(I) » LL I TT( I ) , ST ARN ( MN ) , LLI TT ( MN) , F , Z , RSQAV E ) 

PRINT 104 

CALL RACAH ( LL ITT ( I ) , CAPL ( I ) , LL ITT ( MN ) , C APL ( I ) , CAPLP ( I ) , F , W2) 

PROPER COEFFICIENT CHOSEN FOR TRANSITIONS WITH EITHER 0 CHANGE, 

OR, + OR - 2 CHANGE IN LL IT T. ( SEE EQUATION (C5) ) 

C0EFF=2./3.*(2.*J(N, I ) + 1. ) * ( 2.+CAPLC I )+l . ) **2*RSQAVE 
IF( ABS(LLITT(MN)-LLI TT ( I) ).LT..5> COEFF=COEFF*B IGSML* ( BI GSML+1. ) * < 2 
l.*BIGSML+l.)/( <2.*BI GSML— 1. ) *< 2. *B IGSML+3. ) ) 

IF* ABS (ABS( LLITT*MN)-LLITT( I) )-2. ) .LT. .5)C0EFF=C0EFF*3.*( BIGSML-1. 
1 )*BIGSML/(2.*(2.*BIGSML— 1.)) 

NO=NOJS( MN) 

DO 16 M=1 , NO 

IF( (TJ(M,MN)-TJ(N,I) ).LT.0.)PRINT 107, J(M,MN) 

IF( (TJ(M,MN)-TJ(N,I) ) . L T.O. ) GO TO 16 
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C FURTHER QUADRUPOLE SELECTION RULES APPLIED IN THE NEXT 5 STATEMENTS 

C 

I F( ABS ( J ( N* 1 1 — J ( Mt MN I J • GT *2 • 5 ) GO TO 16 
I F ( ABS ( J ( N* 1 1 + J ( Mt MN ))«LT«1,2) GO TO 16 
IF( ABS(NLITT( n-NLITT(MNI) .GT*.5) GO TO 751 
IF< ABS(LLITT< I)-LLITT(MN) J .GT.,5) GO TO 751 
IF(ABS(J(NfIl-J(MtMN)l*LT#*41 GO TO 17 
751 CALL RACAH (CARL (1) « J<N«n»CAPL(I) « J< M* MN) * CAPS ( 1 1 «F«W1 1 
C 

C QUADRUPOLE TRANSITION PARAMETERS CALCULATED AND PRINTED OUT 

C 

S=COEFF*Wl**2*W2**2* (2* ( M * MN) +1 * ) 

IF(TJ(N, H.EQ.T J(M, MN) I PRINT 115 
IF( TJ(N, n-EQ-TJ(M,MN) ) GO TO 222 
WAADP=2,*PI*C*(TJ(Nt I)-TJ(M f MN) I *<-!•! 
XLAM=l,0E8MTJ(NtI)-TJ(MtMN))*(-l,) 

PRINT 105t J(M t MN) tS, WAADPtXLAM 
I F( OPFV A , NE • 3 ) GO TO 16 

AAAPN=ALPHA*( AOBOHR/C) **4/(320, *( 2,*J(MtMN)+l. ) ) *W AADP**5*S 
PRINT 126 t A A APN 
GO TO 16 

222 PRINT 119 * J ( M* MN) t S 

16 CONTINUE 

17 CONTINUE 
C 

c 

C DELTA CAPITAL L = 1 QUADRUPOLE TERMS 

C 

C 

PRINT 106 
DO 310 M= 1 * NOLEV 

IF( ABS(NLITT( I ) -NLI T T < M ) ) .GT. NL ITO > GO TO 310 
C 

C QUADRUPOLE SELECTION RULES APPLIED IN NEXT 8 STATEMENTS 

C 

IF< ABS(PARITY( I )-PARITY(M)) .GT..5) GO TO 310 
IF< ABS(CAPS( I)-CAPS< M) ) .GT. . 5 ) GO TO 310 
IF(ABS<ABS(CAPL(I)-CAPL<Mn-l.)*GT.*5) GO TO 310 
IF< ABS(LLITT< I ) +LLI T T( M ) ) .L T. . 5 ) GO TO 310 
IF( ABSCCAPSPID-CAPSPIM)) .GT..5) GO TO 310 
IF( ABS(CAPLP( I )-CAPLP<M)).GT..5) GO TO 310 
IF< ABS(LLITT(I)-LLITT(M)).GT.2.5) GO TO 310 

IF(0.5.LT.ABS(LLITT( 1 1 - LL I TT ( MH . AND. ABS ( LL I TT (I) -LL I TT ( M) ) .LT.1.5 
1) GO TO 310 

C 

C FINAL QUADRUPOLE LEVEL IDENTIFICATION PRINTED OUT 

C 

PRINT 110 

PRINT 102 1 NLITT(M»,LLITT(M),CAPS(M> t CAPL<M),PARITY(M) 

PRINT lilt CAPSP(M) tCAPLP(M) 

C 

C DISTINCTION MADE BETWEEN LEVELS WITH SMALLER AND LARGER LLITT 

C 

B IGSML=LL I TT ( I I 

IF( LLITT ( I)-LLITT(M) .LT.O.) BI GSML=L L ITT ( M ) 

IF(ABS(BIGSML-LLITT( III .LT..5J LGR- I 
I F( ABS < BI GSML-LL I TT ( M) ) .LT..5) LGR= M 
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I F ( ABS< BIGSML-LL ITT ( I) I .LT..5I LSM=M 
I F( ABSIBIGSML-LLITTI M) ) .LT..5) LSK=I 
DIFF=STARNILSM)— STARNILGR) 

IF< NPROP.NE. I) GO TO 904 


COULOMB APPROXIMATION QUADRUPOLE RADIAL MATRIX ELEMENT TABLE 
ENTRY POINTS PRINTED OUT 

PRINT 113, ST ARN ( LGR I 
PRINT 114, ST ARN( LSM I 
PRINT 125, DIFF 

SQUARE OF QUADRUPOLE RADIAL TRANSITION INTEGRAL CALAULATED 

904 CALL BATESISTARNI I) , LLITTI I ) .STARNIMI ,LLITTIM1 ,F, Z,RSQAVE> 

PRINT 104 

CALL RACAH l LL ITT 1 1) , CAPL I I ) , LL I TT IM ) , CAPL I M) , C APLP II) , F , W3 ) 

PROPER COEFFICIENT CHOSEN FOR TRANSITIONS WITH EITHER 0 CHANGE, 

OR, + OR - 2 CHANGE IN LLITT. ( SEE EQUATION (C5) ) 

C0EFF=2./3.*I2.*J(N, I> + 1 . ) * (2 . *C AP L (I) +1 . ) *( 2. *CAPL ( M J+ 1. I *W3**2*R 
1SQAVE 

I F ( ABS(LLITT(M)-LLITT(I )) . L T. . 5 ) CO EFF=COEFF*B I GSML* I BI GSML+ 1. ) *12 
1 • *B IGSML+1. I / I I2.*BIGSML-l.l*(2.*BIGSML+3.) J 
I F ( ABS(ABS(LLITT(M)-LLITT( I ) J -2- I . LT . . 5) COEF F=C0EFF* 3. * ( B IG SML-1 . ) 
1 *BIGSML/(2.*12.*BIGSML-1.) ) 

NO=NO J S ( M I 
DO 33 NX= 1 , NO 

IF( (TJ( NX, MJ-TJ(N,I ) J.LT.O. ) PRINT 107, JINX, M) 

IF! (TJINX,M)-TJ(N,I) I.LT.O. )G0 TO 33 


FURTHER QUADRUPOLE SELECTION RULES APPLIED 


IN THE NEXT 5 STATEMENTS 


752 


I F ( ABSI J(NX,M)-J(N, I 
IF! ABSI JINX, MJ+JIN, I 
IF! ABSINLITTI II-NLIT 
IF! ABSILLITTI I )-LLIT 
IFIABSIJI N, I )-J(NX,M 
CALL RACAH (CAPLI I) , 


)).GT.2.5) GO TO 33 
)).LT.1.2» GO TO 33 
TIM IJ.GT..5) GO TO 752 
TIM ) I.GT..5I GO TO 752 
) ) .LT..4J GO TO 310 

JIN, I I, CAPL I Ml ,J(NX,M) , CAPS! I ),F,W4) 


QUADRUPOLE TRANSITION PARAMETERS CALCULATED AND PRINTED OUT 


S=C0EFF*W4**2*(2.*Jl NX,M)+1.) 

IF(TJIN,I)*EQ.TJ(NX,M))PRINT 115 
IFITJIN, II.EQ.TJINX, M) ) GO TO 333 
WAADP=2.#PI*C*(TJIN, I >- TJ I NX , M ) I* ( -1 . » 

XLAM=1.0E8/(TJ( N, I » — T J I NX , M J ) *t -1 . I 
PRINT 105, JINX, M) ,S, WAADP.XLAM 
IFtOPFVA.NE.3I GO TO 33 

A AAPN= ALPHA* I AOBOHR/C) **4/ 1 320. *t 2.* J I NX , M 1+1 . I) *WAADP**5*S 
PRINT 126 , AAAPN 
GO TO 33 

333 PRINT 119, JINX, M) ,S 
33 CONTINUE 
310 CONTINUE 
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C 

C DELTA CAPITAL L = 2 QUADRUPOLE TERMS 

C 

C 

PRINT 108 
DO 31 M=l, NOLEV 

IF( ABS(NLITT( 1 ) -NLI T T( M 1) . GT. NL ITO I GO TO 31 
C 

C QUADRUPOLE SELECTION RULES APPLIED IN NEXT 8 STATEMENTS 

C 

IF( ABS(PARITYd)- P ARI TY( M) ) .GT. . 5 I GO TO 31 
IF! ABS(CAPS(I1-CAPS< M) 1.GT..5) GO TO 31 
IF(ABS(ABS(CAPL( Il-CAPL (Mll-2. 1 .GT..5 1 GO TO 31 
IF(ABS(LLITT(I)+LLITT(M)).LT..5I GO TO 31 
IF( ABS(CAPSP( I ) —CAP S P( M ) ) .GT..5) GO TO 31 
IF( ABS(CAPLP( I )-CAPLP(M)J .GT..5) GO TO 31 
IF( ABS(LLITT( Il-LLI T T( M ) ) .GT. 2. 5) GOTO 31 

IF(0.5.LT. ABS(LLITT( I ) -LL ITT ( Ml) . AND. ABS ( LL ITT ( I) -LL I TT ( M) 1 .LT.1.5 
1) GO TO 31 

C 

C FINAL OUADRUPOLE LEVEL IDENTIFICATION PRINTED OUT 

C 

PRINT 110 

PRINT 10 2 » NLITT(M) , LL I TT < M ) ,CAPS ( M 1 ,CAPL < M 1 , PARI TY ( Ml 
PRINT 111, CAPSP(M) , CAPLP(M) 

C 

C DISTINCTION MADE BETWEEN LEVELS WITH SMALLER ANO LARGER LLITT 

C 

BIGSML=LLITT( I 1 

I Ft LLITTt I) -LLITT (Ml .LT.O.) BIGSML=LL ITT ( M 1 
IFIABSIBI GSML -LLITT (1)1 .LT..5) LGR= I 
IF( ABS(BIGSML-LLITT< Ml ) .LT..5) LGR= M 
I F C ABS(BIGSML-LLITT( I) ) .LT..5) LSM=M 
I Ft ABS(BIGSML-LLITT( Ml ) .LT..5) LSM=I 
DIF F*STARN( LSM1-STARNI LGR) 

I Ft NPROP.NE.l) GO TO 903 
C 

C COULOMB APPROXIMATION QUADRUPOLE RADIAL MATRIX ELEMENT TABLE 

C ENTRY POINTS PRINTED OUT 

C 

PRINT 113, STARN(LGR) 

PRINT 114, STARN(LSM) 

PRINT 125, DIFF 

C 

C SQUARE OF QUADRUPOLE RADIAL TRANSITION INTEGRAL CALAULATED 

C 

903 CALL BATES(STARN( I) , LLITT(I) ,STARN(M) ,LLITT(M) ,F,Z,RSQAVE) 

PRINT 104 

CALL RACAH ( LL ITT ( I 1 , C A PL ( I 1 , LL ITT (M 1 , CAPL ( M ) , CAPLP ( 1 1 , F , W5 1 

C 

C PROPER COEFFICIENT CHOSEN FOR TRANSITIONS WITH EITHER 0 CHANGE, 

C OR, + OR - 2 CHANGE IN LLITT. ( SEE EQUATION <C5) ) 

C 

C0EFF=2.*(2.*J(N,I)+1.)*(2.*CAPL( I )+ 1. ) * ( 2. *CAPL ( M 1 + 1. ) *W5**2*R SQA 
1VE/3. 

I F t ABS(LLITT(M)-LLITT(I 11 .L T. . 5 ) COEFF=COEFF*BI GSML*( BI GSML+1. 1 *(2. 

1 *B IGSML+1 • 1 / ( (2.+BIGSML-1. ) *( 2.*B IGSML+3. ) ) 
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IF(ABS(ABS(LLITT(Ml-LLITT(Ill-2.1 .LT. . 5 IC0EFF=C0EFF*3. * I BIGSML-1 . 1 
I *BIGSML/(2.*(2.*BIGSML-l.l * 

NO® NOJ St Ml 
DO 699 NX=1 , NO 

I F ( (TJ(NX,MI-TJ(N,I1 1.LT.0. 1 PRINT 107, J< NX, Ml 
IFI (TJINX,Ml-TJ(N,Il 1.LT.0. 1G0 TO 699 

FURTHER QUADRUPOLE SELECTION RULES APPLIED IN THE NEXT 5 STATEMENTS 

IF( ABS ( J(NX,MI— J(N, I ) ) .GT.2.51 GO TO 699 
I F( ABS (J(NX,M1+J(N, III.LT.1.21 GO TO 699 
IFI ABSINLITTI Il-NLITTIM I1.GT..51 GO TO 753 
IFIABS(LLITT( I 1—LLI TTI M 11.GT..51 GO TO 753 
IF(ABS(J(N,I1-J(NX,M11 .LT..41 GO TO 31 
753 CALL RACAH I C APL I II , J( N , 1 1 , C APL I Ml , J I NX ,M1 , CAPS! 1 1 , F , W6 1 

QUADRUPOLE TRANSITION PARAMETERS CALCULATED AND PRINTED OUT 

S=C0EFF*W6**2*(2.*J( NX,M1+1.1 
IF(TJ(N» 1 1 . EO. TJ I NX, MI 1 PRINT 115 
IF(TJ(N,I).EQ.TJ(NX,M»IGO TO 444 
WAADP=2.*PI*C*(TJ(N, I)-TJ(NX,M) 1*1-1. 1 
XLAM=1.0E8/(TJ(N,I1-TJ(NX,M1 1*1-1. I 
PRINT 105, JINX, MI, S,WA A DP, X LAM 
I F( OPFVA. NE.3 I GO TO 699 

AAAPN= ALPHA* ( AOBOHR/ C 1**4/ C320.*(2.*J(NX,M)+l.))*WAADP**5*S 
PRINT 126, AAAPN 
GO TO 699 

444 PRINT 119, JINX, MI, S 
699 CONTINUE 
31 CONTINUE 


DIPOLE TERMS 


760 CONTINUE 
F = l. 

PRINT 112 
DO 50 M=l, NOLEV 

IFIABSINL ITT (II— NLITTIMII.GT.NLITOI GO TO 50 

DIPOLE SELECTION RULES APPLIED IN THE NEXT 7 STATEMENTS 

IFI ABSICAPLPI 1 1 — CAPLPI Mil .GT..5I GO TO 50 
I F ( ABSICAPSPI I I-CAPSPIMI I .GT..5I C-0 TO 50 

IFCABSCABSIPARITYIII —PARI TY (MI 1—2. I.GT..51 GO TO 50 
IFIABSICAPLC II-CAPLI Mil .GT. 1.51 GO TO 50 
IFIABSICAPLIII+CAPLIMII .LT.0.5I GO TO 50 
IFI ABS ICAPSIII -CAPS (Mil .GT..5I GO TO 50 
IF(ABS(ABS(LLITT(II-LLITT(MI)-1.I .GT..51 GO TO 50 

DISTINCTION MADE BETWEEN LEVELS WITH SMALLER AND LARGER LLITT 

BIGSML=LLITT(II 

IFI LL ITT! II-LLITTIMI .LT .0. I BIGSML=LL I TT ( MI 
IFI ABSIBIGSML-LLITTI III -LT..51 LGR= I 
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C 

C 

c 


c 

c 

c 

c 


IF( ABS( BIGSML-LLITT (Mil .LT..5I LGR= M 
I F< A8S ( BI GSML-LL ITT ( III .LT..5I LSM=M 
IF( ABS(BIGSML-LLITT< Mil .LT..5I LSM=I 

P=9.*STARN(LGR1**2*ABS( STARN(LGR 1**2-BIGSML**2 I /( Z**2*4. *( 4. *BI GSM 
1L**2-1.1 1 

DIFF=STARN(LSMI-STARN( LGRI 

FINAL DIPOLE LEVEL IDENTIFICATION PRINTED OUT 
PRINT 110 

PRINT 102» NLI TT( Ml t LL I TT(M ) »CAPS ( Ml »CAPL( Ml t PARITY (M) 

PRINT lilt CAPSP(M) tCAPLPIMI 
IFINPROP.NE.il GO TO 707 

COULOMB APPROXIMATION DIPOLE RADIAL MATRIX ELEMENT TABLE ENTRY 
POINTS PRINTED OUT 

PRINT 113 1 STARN(LGR) 

PRINT 114, STARN(LSM) 

PRINT 125tDIFF 
707 $TARNL=ST ARN(LGR) 


C 

C SQUARE OF DIPOLE RADIAL TRANSITION INTEGRAL CALCULATED 

C 

CALL BATES! STARNI I) « LLITTI 1 1 »STARN(M) ,LLITT (Ml »F, Z.RSQAVE) 
PHIS0=RSQAVE*4.*Z**2/(9.*STARNL**2*A8S(STARNL**2-BIGSML**2) ) 

708 PRINT 116 

CALL RACAH ( LL ITT ( 1 1 , CAPL ( 1 1 , LL ITT (M I ,CAPL ( Ml , CAPLP! 1 1 , F , W7 I 
NO=NO JS ( Ml 
DO 49 NX = 1, NO 

I F ( ( T J ( NX tMI-TJCNtll I.LT.0.1 PR I NT 107,J(NX,M1 
IF( (TJINXtMl— TJINtI I I.LT.O.J GO TO 49 
C 

C FURTHER DIPOLE SELECTION RULES APPLIED IN THE NEXT 2 STATEMENTS 

C 

IF(ABS(J(Nt 1 1 + J { NX. Mil. LT. 0.51 GO TO 49 
I F ( ABS ( J ( N. 1 1 — J ( NXt M 1 1 .GT.1.5J GO TO 49 

CALL RACAH C CAPL ( II , J ( N , 1 1 , C APL ( Ml , J ( NX, MI t CAPS (II , F , W8 ) 

C 

C DIPOLE TRANSITION PARAMETERS CALCULATED AND PRINTED OUT 

C 

46 SRAAP=(2.*J(NX, Ml+l. I *( 2. *C APL ( I I + 1 . I *( 2 .*C APL (M I* 1. I *BI GSML*( 4. *B 
1IGSML**2-1. I*W7**2*W8**2 
PR AAP= SR A AP*P 
RAAP=PHISO*PRAAP 
S AAPN=RAAP *(2.*J(N,II+1.1 
I F ( T J ( N , I I.EQ.TJ(NX,MI IPRINT 115 
IF(TJ(N,I1.EQ.TJ(NX,M1I GO TO 555 
WAAP=2.*PI*C*(TJ(N,I I— TJ(NX,MII*(— 1.1 
XLAM=.1E9/(TJ(N,II-TJ(NX,MI 1*(-1. I 
700 PRINT 105, J(NX,MI, SAAPN, WAAP,XLAM 
I F( OPFVA. EQ. 2 1 GO TO 49 

FAAPN=2.*A0B0HRAC3.* ALPHA*C*(2.*J(N, I 1 + 1.1 l*WAAP*SAAPN 
A AAPN=4.*ALPHA*A0B0HR**2/(3.*C**2*(2.*J( NX, Ml+l. II *WAAP**3* SAAPN 
PRINT 717, F AAPN , AAAPN 
GO TO 49 

555 PRINT 119, J(NX, Ml, SAAPN 
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49 CONTINUE 

50 CONTINUE 
15 CONTINUE 

915 CONTINUE 
PRINT 916 
STOP 
ENO 
C 
C 

C RACAH COEFFICIENTS CALCULATED ( SEE EQUATION (10) ) 

C 

c 

SUBROUTINE R AC AH( A, B ,C t D, E« F , W I 
Id* 0. 

CALL DELT A(A«8vE«Dl) 

CALL DELTA(B,D,F,D2) 

CALL DELTA( A,F,C*D31 
CALL DELTA(EtD,CtD4) 

SUM=0* 0 
DO 10 N=1 ♦ 20 
X=N 

Y=X+* 000001 
T 1= Y— A-B- E 
T 2= Y— B-D-F 
T3=Y- A-F-C 
T 4= Y-E-D-C 
T5= A+B+C+D-X 
T6= A+D+E+F-X 
T7= B+C+E+F-X 
I F( T1 • LT.-* 5 ) GO TO 10 
I F( T2.LT.-.5IG0 TO 10 
I F ( T3.LT.-.51G0 TO 10 
I F f T4.LT.— *5)G0 TO 10 
I F ( T5 • LT . — . 5 ) GO TO 1 1 
IF( T6 .LT 5 ) GO TO 11 
I F ( T7.LT.-.51G0 TO 1 1 

DENOM= G AMMF (Tl+i.*0)* GAMM F ( T2+ 1 . ,0 ) * GAMMF( T3+1 . * 0 ) * 

1 GAMMF ( T4+1. » 0)* GAMMF ( T5+1 • ♦ 0 ) * G AMMF ( T6+ 1. , 0) *GAMMF ( T7+1. 
XNUM= (-1. )**N*GAMMF (X+2. «0I 
SUM=SUM+XNUM/DENOM 

10 CONTINUE 

11 K=A+B+C+D 

W=(-l. )**K*D1*D2*D3*D4*SUM 

RETURN 

END 


0 ) 


THE DELTA FUNCTION CALCULATED ( SEE EQUATION (11) ) 


SUBROUTINE DELTA( EM, EN, EL ,D ) 

D=0. 

I F( EM+EL— EN. LT.— .5)G0 TO 20 
I F( EN+EL— EM. LT • — . 5) GO TO 20 
IF(EM+EN-EL.LT.-.5)G0 TO 20 

GT= GAMMFC EM+EN-EL+1 . t 0 )*GAMMF( EN+EL-EM+1* t 01* GAMMF < EL+ EM- EN+1 . 
1 tO) /GAMMF( EM+EN+EL+2.,0) 
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I F ( GT • LT* 0* } WR I TE( 6 * 10 I EM,EN,EL 
10 FOR MAT ( I8H IN DELTA GT NEG,3E17.8) 
D= SORT ( GT I 
20 CONTINUE 
RETURN 
END 


SOUARE OF THE COULOMB APPROXIMATION RADIAL TRANSITION INTEGRAL 
CALCULATED ( SEE EQUATIONS (61,(8), (9) ) 


SUBROUTINE B AT ES ( ST A R , L , STARP , LP , ALF A , Z , RSQU AR » 

DIMENSION A ( 50 ) , AP( 5 0) 

REAL L , LP 

I F ( ST AR-L • LT • • 2 ) GO TO 15 
IF( STARP-LP.LT..2) GO TO 15 
MAX=STAR+STARP+ALFA-.99 
S = S T AR+ST ARP 
A( 11=1. 

AP ( 1 ) = 1 • 

I F ( MAX. LT • 2 ) GO TO 11 
DO 10 1=2, MAX 
XI=I-1 

A ( I )=A(I-1)*(-1.)*STAR*(STAR-L-XI)*( STAR+L-XI+ 1. ) / ( 2 .*X I ) 

APC II=AP( I-ll*(-l.) *STARP*< STARP-LP-X I ) * < ST ARP+LP-X 1+ 1„ ) / ( 2 .*X I) 

10 CONTINUE 

30 GT= GAMMF(STAR+L+1.,0)*GAMMF(STAR-L,0)*GAMMF(STARP+LP+1.,0) 

1 *GAMMF( $TARP-LP,0) 

I F ( GT*LT# 0* ) WR ITE ( 6 ,20 ) STAR, L, STARP 
20 FORMAT ( H IN BATES GT NEG*3E17.8J 

11 F=2.**S*STAR**( STARP+ALFA+2. )*ST ARP** < STAR+ALFA+2. >/($**( S+ALFA+3 • 
1)*Z**ALFA* SQRT(GT) ) 

MAX=MAX+1 

14 $UM=0. 

DO 12 K= 1 ♦ MAX 

DO 13 KP= 1 , MAX 

IF( K+KP.GT.MAX) GO TO 12 

I F ( STAR+STARP+3.+ALFA-FL0AT(K+KP).LT..2) GO TO 15 

SUM=SUM+A(K)*AP(KP) *(1. /STA R+ 1. /ST ARP ) ** ( K + KP ) *GAMMF ( S+'ALF A + 3.-FL0 
1 AT( K+KP ) , 0) 

13 CONTINUE 

12 CONTINUE 
SUM=SUM*F 
RSQUAR=SUM**2 
GO TO 16 

15 RSQ AR=0. 

16 RETURN 
END 
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Sample Input Deck 


CARO 1 ( OPTIONS I 

OP QUA, QPFVA j NPROP , NL ITO 

2 3 1 2.0 


CARD 2 ( ATOMIC SYSTEM DEFINED, HELIUM I I 

EL EMNO,Z, NOLEV, NOSER 

2 1. 6 I 


CARD 3 ( THERE IS ONLY 1 CARD LIKE CARD 3, SINCE HE I HAS ONLY 1 SERIES ) 

CAPSSI I) ,CAPLSm , SLIMII) 

.5 0. 198305.7 


THE REMAINING DATA CARDS ARE IN TWO-CARD SEQUENCES WITH ONE TWO-CARD 
SEQUENCE FOR EACH INPUT LEVEL I THUS THERE ARE SIX TWO-CARD SEQUENCES, 
SINCE NOLEV = 6 I 


1ST CARD OF 1ST SEQUENCE 

NLITT ( II ,LL ITT (11, CAPS (II, CAPL(l), PARITY (II, CAPSPIll ,CAPLP(1I 

2. 0. 1. 0. 2. 0.5 0. 

2N0 CARD OF 1ST SEQUENCE 
T J (N ,1 l , N=1 j NO J (1ST LEVEL I 

159850. 318 


1ST CARD OF 2ND SEQUENCE 

NL ITT( 21 ,LLITT (2 I , CAPS(2I,CAPL (2 I, PARITY (21 ,-CAPSP ( 21 ,CAPLP(2I 

2. 1. 1. 1. 0. 0.5 0. 

2N0 CARD OF 2ND SEQUENCE 
T J ( N , 2 I , N =1 , NO J (2ND LEVEL I 

169082.185 169081.189 169081.111 
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C 1ST CARO OF 3RD SEQUENCE 

C NL ITT ( 3» ,LLITT (31 , CAPS ( 3 > , CARL ( 3 » , PARITY ( 3 ) , CAPSP ( 3 ) , CAPLP ( 3» 

3. 1. 1. 1. 0. 0.5 0. 

C 2ND CARD OF 3RD SEQUENCE 

C TJ (N,3I ,N=1,N0J (3RD LEVEL ) 

185559.277 185559.085 185558.920 

C 1ST CARD OF 4TH SEQUENCE 

C NLITT(4t ,LLITT<4> , CAPS (A* , CARL (4) » PARITY (4), CAPSP (4> , CAPLP (4> 

3. 2. 1. 2. 2. 0.5 0. 

C 2ND CARO OF 4TH SEQUENCE 

C TJ (N ; 4» ,N=1,N0J <4TH LEVEL * 

186095.900 186095.900 186095.900 

C 1ST CARD OF 5TH SEQUENCE 

C NL ITT ( 5 ) ,LLITT(5» , CAPSC5) , CAPL ( 5 > , P AR I TY ( 5 ) , CAPSP (5 » , CAPLP (5 ) 

4. 1. 1. 1. 0. 0.5 0. 

C 2ND CARO OF 5TH SEQUENCE 

C TJ <N r 5) ,N=1,N0J C5TH LEVEL ) 

191211.420 191211.420 191211.400 

C 1ST CARD OF 6TH SEQUENCE 

C NLITT(6)jLLITT(6)>CAPS(6)jCAPL(6)»PARITY(6)»CAPSP(6) .CAPLP (6* 

4. 2. 1. 2. 2. 0.5 0. 

C 2ND CARD OF 6TH SEQUENCE 

C TJIN,61,N=L,N0J L6TH LEVEL ) 

191438.830 191438.830 191438.830 
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Sample Printout 


c 

C HELIUM I 

C OPTION VARIABLE VALUES ARE OPQUA=2 ,0PFVA*3 y NPROP=l y NL IT0=2 

C 


ELEMENT NUMBER t 
INITIAL LEVEL 
LITTLE N= 2* 
PARENT TERM 
INITIAL J VALUE 


2, WITH Z= 

LITTLE L= 0* 
CAPITAL S= 

J= l.Oy 


1 

CAPITAL S= 1*0, CAPITAL 
• 5 y CAPITAL L* 0.0 
RAA* 21.790 


L= 


0.0 


DELTA L EQUALS 0 QUADRUPOLE TERMS 


DELTA L EQUALS 1 QUADRUPOLE TERMS 


DELTA L EQUALS 2 QUADRUPOLE TERMS 


FINAL LEVEL 


LITTLE N= 3, 

LITTLE L= 2, CAPITAL S= 1.0, CAPITAL 

L= 2.0 

PARENT TERM 

CAPITAL S= .5, 

CAPITAL L* 0.0 


N 

STAR OF LARGER 

LITTLE L VALUE* 

2.998 


N 

STAR OF SMALLER LITTLE L VALUE* 

1.689 


DIFFERENCE* -1. 

309 



FINAL J VALUES 

S 

WAADP 

LAMDA 


1.0 

2.63 06E+02 

4.9438E+15 

3810.2 

A* 

5.73 19E+00 





2.0 

4.3842E + 02 

4.9438E+15 

3810.2 

A* 

5.73 19E+00 





3.0 

6. 1379E+02 

4.9438E+15 

3810.2 

A* 

5.73 19E+00 





FINAL LEVEL 

LITTLE N* 4t LITTLE L= 2y CAPITAL S= l.Oy CAPITAL L= 2.0 
PARENT TERM CAPITAL S* .5t CAPITAL L* 0.0 

N STAR OF LARGER LITTLE L VALUE= 3.998 
N STAR OF SMALLER LITTLE L VALUE* 1.689 
DIFFERENCE* -2.308 


.0 


t PARITY* 2 


y PARITY* 2 


y PARITY* 2 
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FINAL J VALUES 

s 

WAADP 

LAMDA 

1.0 

1.9120E+00 

3.47 44E+01 

5.9502E+15 

3165*7 

2.0 

= 1.9120E+00 

5.7907E+01 

5.9502E+1 5 

3165.7 

3.0 

= 1.9120E+00 

8. 1069E+01 

5.9502E+15 

3165.7 


DIPOLE TERMS 


FINAL LEVEL 

LITTLE N= 2, LITTLE L= I» CAPITAL S = 1.0, CAPITAL L= 1.0, PARI TY= 0 
PARENT TERM CAPITAL S = .5, CAPITAL L = 0.0 

N STAR OF LARGER LITTLE L VALUE* 1.938 
N STAR OF SMALLER LITTLE L VALUE* 1.689 
DIFFERENCE* -.249 

FINAL J VALUES S WAAP LAMDA 

.0 6.4877E+00 1.7390E+15 10832.0 

ABS. OSCILLATOR STRENGTH* 6.0643E-02 A* 1.0342E+07 


1.0 

ABS. OSCILLATOR 


1.9463E+01 1.7390E+15 

STRENGTH* 1.8191E-01 A* 1.0338E+07 


10832.0 


2.0 

ABS. OSCILLATOR 


3.243 8E+01 1.7390E+15 

STRENGTH* 3.0318E-01 A* 1.0338E+07 


10832.0 


FINAL LEVEL 

LITTLE N= 3, LITTLE L= 1, CAPITAL S= 1.0, CAPITAL L= 
PARENT TERM CAPITAL S= .5, CAPITAL L* 0.0 

N STAR OF LARGER LITTLE L VALUE* 2.934 
N STAR OF SMALLER LITTLE L VALUE* 1.689 
DIFFERENCE* -1.245 

FINAL J VALUES S WAAP 

.0 2. 553 3E-01 4.8426E+15 

ABS. OSCILLATOR STRENGTH* 6.6464E-03 A* 8.7900E+06 


1.0, PARITY* 


LAMDA 

3889.7 


1.0 7.6598E-01 4.8426E+15 

ABS. OSCILLATOR STRENGTH* 1.9939E-02 A* 8.7898E+06 


3889.7 
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2.0 1.2766E+00 4.8426E+15 3889.7 

ABS. OSCILLATOR STRENGTH* 3.3231E-02 A= 8.7896E+06 


FINAL LEVEL 

LITTLE N= 4, LITTLE L= 1, CAPITAL S= 1.0» CAPITAL L= 1.0, PARI TY= 0 
PARENT TERM CAPITAL S= .5, CAPITAL L= 0.0 

N STAR OF LARGER LITTLE L VALUE* 3.933 
N STAR OF SMALLER LITTLE L VALUE= 1.689 
01 FFERENCE= -2.244 

FINAL J VALUES S WAAP LAMDA 

.0 8.2923E— 02 5.9073E+15 3188.7 

ABS. OSCILLATOR STRENGTH* 2.6331E-03 A* 5.1819E+06 


1.0 2. 48 7 7 E- 01 5.9073E+15 3188.7 

ABS. OSCILLATOR STRENGTH* 7.8994E-03 A* 5.1819E+06 


2.0 4. 146 IE- 01 5.9073E + 15 3188.7 

ABS. OSCILLATOR STRENGTH* 1.3166E-02 A* 5.1819E+06 
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RELATIVE LINE AND MULTIPLET STRENGTHS FOR DIPOLE TRANSITIONS 


In the notation of Shore and Menzel (ref. 20), the line strength may be factored as 


s (0,0/) = S(a',a) = R^R*^ I 2 


(Bl) 


where Rj ine and R mult are the line factor and multiplet factor, respectively, and I 
the radial factor. They are related to the customary relative line and multiplet strengths 
(tabulated by Goldberg in refs. 21 and 22) and the radial factor a by 

s(line)= ^T (B2) 


S(mult) = Z>(2S + l)(ul - l)Rmult 


(B3) 


o2 = 



(B4) 


where S is the total atomic spin quantum number and l > is the larger of the individual 
electron orbital angular momentum quantum numbers l and l\ 

The radial factor I is related to the radial transition integral I s by 



Under the restrictions of LS coupling and transitions involving no equivalent elec- 
trons, the relative line and multiplet strengths can be expressed in terms of Racah coef- 
ficients W(abcd;ef) (ref. 9) as 


and 


S(line) = 


(2J + 1)(2J’ + 1) w 2( LJL ,j, s1 ) 

(2S + 1) 


(B6) 


S(mult) = (2S + 1)(2L + 1)(2L' + 1)Z>(4Z 2 - l)w 2 ^LZ'L’;Lil) (B7) 

The usual notation has been used for the total atomic orbital angular momentum quantum 
number L and the total atomic angular momentum quantum number J; the subscript 1 
indicates that the quantum number refers to the parent configuration. Racah coefficients 
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can be found in tabular form (ref. 10) or can be calculated by computer (eq. (10)). In 
order to generalize the results to equivalent electron configurations, coefficients of 
fractional parentage are required (ref. 20). 
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QUADRUPOLE RELATIONS 

For quadrupole transitions (ref. 23), the transition probability can be expressed in 
terms of the quadrupole line strength as 


A(oi,Q!') = 


1 /e 2 \ 6 c 1 ( l T a " Ta’I \ S{a,a') 


320\Kc/ aQ g^t \ Roo 


e2a 0 


(Cl) 


The quadrupole line strength can be factored analogously to the dipole line strength 


as 


R Une R muU l2 < C2 > 

where the angular factors can be expressed as 

R line = (-1) S " L " J ’[(2J + 1)(2J* + 1)] 1//2 W(LJL’J*;S2) (C3) 


and 

R mult = (-l) Ll " i " L ’[(2L + 1)(2L* + ljj^W^LZ’L'jL^j 


(C4) 


The radial factor is given by 

r 

fi (; + 1)1,21 + i ) 




" 1 1/2 


(2 1 - 1)(2 1 + 3) 


12 


3(z< + !)(;< + 2) 
2(2i< + 3) 


n 1 / 2 


h 


(A l = 0 transitions) 


(A l = ±2 transitions) 


where I 2 is the quadrupole radial transition integral I s (defined by eq. (6)) with 
s = 2. 


(C5) 
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SCALED THOMAS-FERMI METHOD 

The scaled Thomas -Ferm_ method of Stewart and Rotenberg (ref. 4) makes use of 
the Thomas-Fermi statistical model of atoms and ions. The potential due to the core of 
a multiple electron atom can be approximated by the potential of an appropriate Thomas- 
Fermi ion at all radial distances (since outside the Thomas-Fermi ion the potential 
becomes the correct Coulombic one), namely a Thomas-Fermi ion with net charge and 
nuclear charge equal to the core charge and nuclear charge of the actual atom or ion, 
respectively. Stewart and Rotenberg introduce a scaling factor which allows for uniform 
contraction and dilatation of the Thomas-Fermi core in an effort to compensate for the 
inadequacies of the Thomas-Fermi model. The scaling factor is determined by requiring 
that the energy parameter be identical to the measured energy value of the system. 
Stewart and Rotenberg have shown that the wave functions obtained by the scaled Thomas- 
Fermi method compare favorably with the wave functions obtained by the more sophisti- 
cated self-consistent-field method. In the region outside the Thomas-Fermi core, the 
scaled Thomas-Fermi method yields radial wave functions which are identical to the 
Bates -Damgaard Coulomb approximation wave functions except for a difference in the 
normalization factor. 

A problem sometimes encountered in applying the scaled Thomas-Fermi method 
deserves mention. A Thomas-Fermi ion exists only when the binding energy of the sys- 
tem is greater than the corresponding hydrogenic value. When this requirement is not 
met, a negative Stewart -Rotenberg scaling factor results which implies that the Thomas- 
Fermi ion has a radius smaller than zero, that is, no Thomas-Fermi ion exists. In such 
cases the Bates -Damgaard Coulomb approximation should be used. 

Such a situation exists in neutral helium for all the singlet states given in table I. 
Instead of abandoning the scaled Thomas-Fermi method for the singlet states, the binding 
energies of the singlet states were replaced by the average energies obtained from 

^ _ E (singlet) + 3E (triplet) /rkiN 

h*av ^ \Di) 

This artifact makes the Stewart-Rotenberg scaling factor positive and finite for all the 
singlet states. The results of these calculations for the singlet states are presented in 
table I (enclosed in brackets), but in general they may be expected to be less accurate 
than the corresponding Coulomb approximation results. 
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